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Abstract 

Recently, a homogeneous superfluid state with a single gapless Fermi surface 
was predicted to be the ground state of an ultracold Fermi gas with spin population 
imbalance in the regime of molecular Bose-Einstein condensation. We study vortices 
in this novel state using a symmetry-based effective field theory, which captures the 
low-energy physics of gapless fermions and superfluid phase fluctuations. This theory 
is applicable to all spin-imbalanced ultracold Fermi gases in the superfluid regime, 
regardless of whether the original fermion pairing interaction is weak or strong. We 
find a remarkable, unconventional form of the interaction between vortices. The 
presence of gapless fermions gives rise to a spatially oscillating potential, akin to 
the RKKY indirect-exchange interaction in non- magnetic metals. We compare the 
parameters of the effective theory to the experimentally measurable quantities and 
further discuss the conditions for the verification of the predicted new feature. Our 
study opens up an interesting question as to the nature of the vortex lattice resulting 
from the competition between the usual repulsive logarithmic (2D Coulomb) and 
predominantly attractive fermion-induced interactions. 
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1 Introduction 



The appearance of quantized vortices is a hallmark of superfluid flow. Vortices 
have been studied for decades, experimentally and/or theoretically, in a vari- 
ety of systems as diverse as type-II superconductors, superfluid Helium liquids, 
rotating ultracold atomic Bose and Fermi gases, and neutron stars. [1] Among 
these systems, the quantum gas of resonantly interacting fermionic atoms with 
equal populations of both (hyperfine) spin components, a prototype system 
for the interesting BEC-BCS crossover physics, [2|3lHll5]|6] has been inten- 
sively studied over the past several years. The first experimental observation 
of vortices in 2005 by the MIT group provided a definitive evidence for super- 
fluidity in atomic Fermi gases. [7J In addition, several theoretical studies have 
analyzed the possible new properties of vortices across a Feshbach resonance 
from the BCS to the BEC side. IgMMTTl 



The physics of atomic Fermi gases is also of fundamental interest beyond 
the standard BCS/BEC physics, owing to the new tuning flexibility in the 
atomic gas systems. Under the condition of density imbalance (hence mis- 
matched Fermi surfaces) between the spin-up and -down fermions, a mod- 
ulated Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) [12] superfluid phase, has 
long been theoretically anticipated. The theoretical interest in pairing with 
mismatched Fermi surfaces has been revitalized by the proposal of breached- 
pairing superfluidity, with a number of exotic superfluid states being proposed 
or revisited. P^pi] Breached-pairing superfluid phase with two gapless Fermi 
surfaces (BP2), related to the unstable Sarma phase, [15] was found to be sta- 
ble under the introduction of new effects, such as the mass imbalance and/or 
momentum-dependent pairing interaction. [16j|17j Important developments in 
the subject are recent studies by various groups [T51I19II20] investigating the 
Feshbach-resonant regime of strong interactions. [2T] The first experiments on 
ultracold fermionic gases with spin population imbalance have recently been 
carried out [221123] and thereby brought the subject to the forefront of the cold 
atom physics. The imbalanced Fermi gas is presently the subject of fervent 
research activities. [24,25, 26, 27 JS] 



One of the states commonly found in various theoretical approaches [T8|19II20] 

is a homogeneous superfluid with a single gapless Fermi surface on the molecu- 
lar (BEC) side of the Feshbach resonance (the BEC regime). This state, which 
consists of coexisting molecular superfluid and fully-polarized Fermi gas of the 
majority-spin component, is closely related to the BP2 phase, [TBPT] but dif- 
fers from the latter in the number of gapless Fermi surfaces. We will refer to 
this phase as BPl (breached-pairing state with a single gapless Fermi surface) 
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after Ref. [25] and [28] (also dubbed as "magnetized superfluid" and denoted 
SFm in Ref. [19] and [21]). The BPl phase is predicted to robustly exist in 
a relatively large area in the phase diagram of the spin imbalanced Fermi 
gases (also called "polarized Fermi gases"). Additionally, several theoretical 
works have found the analog of BPl in a trap, the so-called superfiuid-normal- 
mixture phase. [29] It is then of great interest to examine properties of this 
novel superfluid. 

When fermionic excitations are fully gapped, the physics of vortices belongs to 
the universality class of the XY model where the phase of the superfluid order 
parameter plays the dominant role. The vortex sector of this model is described 
by a 2D Coulomb gas of "charges" with a repulsive logarithmic interaction. 
This paradigm is challenged in a fundamental fashion in the BPl state. The 
presence of gapless fermionic quasiparticles is the distinguishing feature of 
this superfluid phase and is expected to have important consequences for its 
physical properties. Since the observation of a vortex lattice is perhaps the 
only unambiguous signature of superfluidity in ultracold fermionic gases, [7] 
it is of interest to examine possible ramiflcations of the presence of gapless 
fermions on the interaction between vortices in this system. 

In this work, we determine the effective interaction potential between the 
vortices in the BPl phase. Because we are solely concerned with the intrin- 
sic effect of the gapless fermions on the interaction between vortices in the 
BPl state, we consider only the homogeneous case, disregarding the effects 
of traps. We exploit the method of effective fleld theory, [30] based explicitly 
on broken continuous symmetry. [31] This method is particularly suitable for 
problems involving strongly-coupled systems in the long-wavelength limit and 
has already proven to be fruitful in treatments of strongly-interacting regimes 
of ultracold atomic gases. [52|33|34|[55] In the present context, the relevant 
degrees of freedom for a low-energy effective fleld theory are the superfluid 
phase fleld and the fleld describing the gapless fermionic excitations. 

We show that the resulting interaction between vortices in the BPl phase is 
not of the pure Coulomb form, but contains an additional fermion- induced con- 
tribution that oscillates on a length scale set by the spin polarization, closely 
resembling the Ruderman-Kittel-Kasuya-Yosida (RKKY) indirect-exchange 
interaction in non-magnetic metals. [36] In order to show that such an unusual 
vortex interaction is perfectly compatible with the BPl phase, we calculate 
the superfluid density from a microscopic model in the parameter regime rele- 
vant for the BPl state. We demonstrate that the superfluid density in BPl is 
positive throughout, which corroborates the dynamical stability of this phase 
and warrants its further investigation. 

The outline of the remainder of this paper is as follows. In Sec. II we introduce 
notation and conventions for fermion quasiparticles and superfluid phase to be 
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used throughout. In Sec. II wc present self-contained field-theoretical deriva- 
tion of the superfluid density, followed by the calculation of this quantity in 
the parameter regime where the BPl phase is realized. Sec. IV starts with 
a low-energy effective field theory for gapless fermions and superfiuid phase 
field, from which we derive an effective theory for phase fiuctuations by in- 
tegrating out the fermionic degrees of freedom. In Sec. V we first derive the 
effective theory for vortices and their effective interaction in momentum space. 
Then we present the calculation of the effective vortex interaction potential 
in real space, accompanied by the discussion of its physical significance. Fi- 
nally, we summarize in Sec. VI. Some mathematical details are relegated to 
the Appendices. 



2 Notation and conventions 

2.1 Gapless branch of fermionic quasiparticles 

The Bogoliubov-quasiparticle energy spectrum of the system containing two 
fermion gases with equal masses (m) and unequal chemical potentials ^ 
III) is given by (in what follows H—l, unless stated otherwise) 



where /i — (//-f -|- //J/2 is the average chemical potential (thermo dynami- 
cally conjugate to the overall atomic number of two species combined) and 
5 = /i| — measures the mismatch between the two chemical potentials (con- 
jugate to the relative density imbalance). Since our treatment concerns the 
spin-polarized homogeneous superfiuid realized deeply on the BEC side of the 
Feshbach resonance, the average chemical potential is assumed to be negative 
{Ijl < 0) in what follows. For definiteness, we hereafter also assume that S > 0. 

For sufficiently large mismatch S (more precisely, for 6/2 > ^/JP~+~A^), the 
lower branch of the above quasiparticle dispersion is gapless, with a single 
effective Fermi surface. Hereafter, for the sake of brevity, we denote it as 



The effective Fermi surface is defined by = for |k| = A;^, where fc^ is the 
radius of the "breached-pairing Fermi ball" in momentum space. The latter 
is controlled by the density imbalance Uf, = — between the two pseudo- 
spin components (here > n|, as a consequence of the fact that 5 > 0), 




(1) 
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which, as imphed by the Luttinger theorem [37] , is equal to the volume of the 
"breached-pairing Fermi ball" in momentum space 

r d^k ^ _kl_ 

J\moM (27r)3 67r2 ■ ^ ^ 

This leads to a simple expression for kf, : 

h=[Qn'in,-n,)]y' . (4) 

An example of gapless dispersion given by Eq. ([2]) is depicted in Fig. ID.li 



2.2 Superfluid phase field and its decomposition 



The superfluid phase field 6'(x, r) represents the phase of the complex Cooper- 
pair amplitude : 

{^PM = \{^M\e'' , (5) 
where ip^ and ipi are the fields describing fermions of opposite spins. In the 
following, we use the standard decomposition of the superfluid phase field 
into a static (classical) contribution 6't,(x) (singular part) and a quantum- 
fluctuating contribution 0(x, r) (regular part): 

e(x,r) = e,(x)+0(x,r) . (6) 

(/)(x, r) is a non-compact (unbounded) field describing "spin-wave" (smooth, 
i.e., non-topological) phase fluctuations ; 6'„(x) is a multi- valued field pertain- 
ing to the topological defects of broken global U{1) symmetry — vortices. 

Superfluid vortices can effectively be treated as classical point-objects in two 
dimensions. In reality, a two-dimensional theory is applicable to rotating su- 
perfluids as long as the rotation frequency is not too high and the flow is 
everywhere confined to a plane perpendicular to the rotation axis. ^38j These 
conditions allow formation of straight vortex lines parallel to the axis of ro- 
tation. The arrangement is then essentially two-dimensional, equivalent to an 
array of point vortices with circulation of the same sign. 

The gradient of the singular part of the phase field can conveniently be ex- 
pressed as [39] 

Ve, = Ko{e, X V) 1 rf^x' G(x,x')p(x') , (7) 

where kq = h/{2m) = irh/m {h restored for the sake of clarity) is the circu- 
lation quantum, p(x) stands for the vortex "charge density", and G'(x, x') = 
G(|x — x'l) is the Green's function of the two-dimensional Laplacian: 

V'G(x,x') = (5(2nx-x') . (8) 
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The vortex charge density is defined as 

p(x) = 27r^iV,5(2)(x_x,) , (9) 

a 

where Na is the winding number, viz. topological charge, of a vortex located 
at position in the xy-plane and is the unit vector perpendicular to this 
plane. Following the standard prescription, [STIHO] the gradient of 6^ can be 
associated to a vortex gauge (vector) field a(x) through a = —VO^, which by 
virtue of Eq. ([7]) becomes 

a= -fi;o(e, X V) y rfVG(x,x')p(x') . (10) 

Consequently, the vortex gauge field obeys the condition 

V X a = -Kop(x)e^ , (11) 

whose momentum-space version (obtained by spatial Fourier transformation) 
reads 

- X aq = Kop(q)e^ , (12) 

where q is a two-dimensional wave-vector (q ■ = 0) of phase fiuctuations 
and p(q) is the Fourier transform of the vortex charge density p(x). 

We now make use of the decomposition aq = a|^ + a^ of the vortex gauge field 
into the longitudinal and transverse components with respect to the momen- 
tum q, respectively, and they are given by 



„ll _ (q-aq)q „x _ (q X aq) x q 

.1 



Because q x all = 0, Eq. (1121) implies that 



p(q)e^ = q X a^ , p(-q)e^ = — q x a^ , (14) 

«o /to 

where the second equation in f[T^ has been obtained from the first one by 
a simple replacement q — > — q. The scalar product of the last two equations 
yields 

P(q)p(-q) = — (q X a^) • (q X a^^) . (15) 

Kg 

Now using the fact that q ■ = q • Si^^ = 0, together with the identity of 
vector algebra 

(A X B) ■ (C X D) = (A • C)(B ■ D) - (A • D)(B ■ C) , (16) 
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it is straightforward to obtain a useful relation 



q2 

P(q)p(-q) = — (ai ■ a^q) . (17) 



For convenience, we henceforth adopt the Coulomb gauge V ■ a = 0. In mo- 
mentum space this reads q-aq = 0, which means that in this gauge the vector 
field aq is purely transverse (a|[ = 0, viz. aq = a^). Consequently, Eq. ( ITTI) 
can be rewritten as 

P(q)p(-q) = — (aq • a_q) , (18) 
a form that will be used in the following sections. 



3 Superfiuid density calculation 



The superfiuid density is a quantity of paramount importance in the realm 
of quantum liquids. This macroscopic observable has the nature of a trans- 
port coefficient and describes the response of a superfiuid system to a Galilean 
boost transformation. Its low-temperature behavior reflects the key properties 
of the ground state. [H] In what follows, we calculate the superfiuid density 
in the BPl state, to be subsequently used as an input to the effective field 
theory in the second part of this paper. In order to make the presentation 
self-contained, we start from a microscopic fermion-pairing model and de- 
rive a general expression for the superfiuid density by following the standard 
field-theoretic method of Ref. [12]. This approach has proved to yield equiva- 
lent results as that of Ref. We then specialize to the case of equal mass 
fermions and evaluate the superfiuid density in the relevant parameter regime 
for the realization of the BPl phase. 

Our starting point is the microscopic Lagrangian 

Co= E C f 5r - ^ - P.) + #T>I^i^T (19) 



describing pairing of two species (denoted by a formal pseudo-spin variable 
=1) i) of fermions with masses rricr, chemical potentials fia, and attractive 
inter-species contact interaction with coupling constant g. In the mean-field 
approximation, the thermodynamic potential for this system is given by 
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where 

G'=\ ' I (21) 

is the inverse of the fermion propagator in the Nambu space, with free fermion 
dispersions ek,o- = k^/ {2m^)—fi„ and A set to be reaL This fermion propagator 
has the matrix form 



(22) 



with 



and e^, defined as 



6^ = - (6k,r ± eka) , eA = + A2 • (27) 

Quasiparticle energy spectrum is determined by the poles of propagator (12^ . 
i.e., by the solution of equation det Q^^ = : 

^k=eA + ek , E^=eA-e^. (28) 



The occupation numbers of two species of fermions can be calculated from the 
diagonal elements of propagator 



nt(k)=r'limE6;ne^""' , (29) 



lUJn 



n^{k) = -p-' hmY^g^^e-'^-^ . (30) 
Upon performing Matsubara frequency summations we obtain 



n^{k)=uinpiE^) +vinp{-E^) , (31) 
n^{k)=uinF{E^)+vlnF{-E+) , (32) 
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where npiz) = {exp{f3z) + 1) ^ is the Fermi distribution function and 

are the coherence factors (squared Bogohubov amphtudes). 

Under Gahlean boost with velocity v^, A transforms as A ^ /\g*(™i+™-T)'*'s-x^ 
while fermion fields transform as ipa- — > ipa^^'^'"'^, where qo- = m^Ws. The 
superfluid (mass) density tensor pij is defined through 

= fi(0) + j, ■ V, + lp,^.(v,),(v,), + 0{yI) . (34) 

For a homogeneous and isotropic superfluid this tensor is diagonal, viz. pij = 
5ijPs/3, where ps is the superfluid mass density. Accordingly, the last formula 
reduces to 

fi(v,) = fi(0) + js ■ V, + ^p,v2 + 0(v3) . (35) 

o 

By transforming the fermion fields and A according to the rules stated above, 
the thermodynamic potential becomes 

l^(v.) = ---/5-E/7^trln^- , (36) 
where Q~^{iujn,^) is the v^-dependent fermion inverse propagator 

A icUn + ek-m|Vs,i 

It is easy to check that the latter can be expressed as 

g;' = g-' - (k ■ v,)i2x2 - ^v^s^ , (38) 

where = diag(m-f, —mi), i.e., as 

q;' = {i - (k- v,)^ - Ul{g^m)] . (39) 

By making use of the well-known expansion formula ln(l— 2;) = z+^z'^+0{z^), 
we obtain 

trlng;;! = trln^-i + (k-v,)tr(e;)-^tr(^S„)-i(k-v,)2tr(6;2)+0(v3). (40) 

The last result, combined with Eq. (l36ll . enables us to expand the thermody- 
namic potential fl^Vg) in powers of and read off the superfluid density from 



9 



the quadratic term 



Ps = m^rii + m^n^ + 



where 



(41) 



(42) 
(43) 
(44) 



By carrying out these Matsubara frequency summations we get 



nF{E^)+nF{E^) 

nF{E^)+nF{E^) 

l-nF{Et)-nF{E-^ 
ca 



^k^k + ^F(^k>k + <(^k)4 



2 2 
«k% , 



(45) 
(46) 
(47) 



where n'p{x) = dnp^x) / dx. Using the last three equations and identity + 
= 1, we then obtain 



d^k k2 r , 



(27r)3 3 



n'piE^) + n'p{E^ 



(4J 



Here and are momentum-space integrals of ^^(k) and ?T-|(k), respectively. 
In what follows, we employ formula fHSl) to determine the superfiuid (number) 
density Ug = Ps/ {m-^ + mi) at zero temperature in the special case of an equal 
mass system of interest in the present work. 

At zero temperature np{x) = 6{—x), whereby for = mj^ = m we readily 
obtain 

d^k k2 r 



Ps = m{n^ + 72°) 



6{E^) + 6{E'^ 



(49) 



(27r)3 3 

where and are the zero-temperature values of and n^, i.e., the respec 
tive momentum-space integrals of 



n',{k)=uld{-Et)+vid{E^\ 
nl{k)=ul9{-E^)+vle{Et\ 



(50) 
(51) 
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The squared Bogoliubov amplitudes (coherence factors) in this special case 
are given by 



After performing a trivial angular integration, on account of the fact that the 
upper branch is always positive in our case and that = £k, we obtain 
an expression for the superfiuid density Ug = Ps/ (2m) : 



We now invoke the property of Dirac's S function 

5(/(x)) = ^^E^(^-^^) ' (55) 

Xi being the simple zeros of the function f{x) (i.e., f{xi) — 0, f'{xi) ^ 0). 
This property can be equivalently stated as 

/^AWWW)<fc = i:^. (56) 

where the last sum extends over all the simple zeros of f{x) within the domain 
of integration D. On account of the fact that |k| = /c^ is the only zero of the 
function £(|k|), simple transformations lead to 



1 P \/(|^-//)2 + A2 

As can straightforwardly be derived, n? and n? are given by 



1 /""^ 



Their calculation requires numerical evaluation of the integral 
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Finally, we eliminate kb from Eq. (1571) (in favor of parameters fi, S, and A) 



The superfluid density is calculated numerically based on the derived expres- 
sions. Some typical results thereby obtained for the superfluid density as a 
function of the pairing gap are presented in Fig. ID.2[ In Fig. ID. 31 superfluid 
density is plotted as a function of the spin-polarization P = (rij^— n|)/(n| + ?T,|) 
for different values of the pairing gap. In contrast to the related BP2 state, no 
anomalous negative value of the superfluid density in the BPl phase is found, 
thus corroborating the dynamical stability of this phase. 



4 Low-energy effective field theory 

In this Section, we start from a symmetry-based low-energy effective La- 
grangian for the gapless branch of fermions and the superfluid phase field. As 
is widely accepted, deep in the superfluid regime the dominant role is played 
by the superfluid phase fluctuations, while the fluctuations of the amplitude 
of the order parameter can be neglected (the London limit). [39] We then de- 
rive an effective phase-only action by integrating out the fermionic degrees of 
freedom. The upper cutoff for the wave vector q of the phase fluctuations is 
set by = {2m A)^^"^, a momentum scale corresponding to the pairing gap 
A. This is a consequence of the BPl superfluid phase being realized in the 
strong-coupling regime on the BEC side of a Feshbach resonance, where the 
pairing gap is related to the binding energy of a Feshbach molecule. In more 
common examples of fermionic pairing (e.g., the weak-coupling regime on the 
atomic side of a Feshbach resonance) the momentum cutoff would have been 
set by ,^0^"'^ ~ A/wi? (the inverse of the coherence length), where vp is the av- 
erage Fermi velocity of the pairing fermions. It is important to point out that 
the magnitude |q| of the wave vector of phase fluctuations is not necessarily 
small as compared to kb- Namely, as already stated in Sec. 12.11 (recall Eq. (jl])), 
kb is controlled by the spin population imbalance and can therefore (by tuning 
the population imbalance) be made arbitrarily small. 

The form of our effective theory will be chosen so as to obey the Galilean 
invariance. In a Galilean invariant system the momentum density Toi (the off- 
diagonal part of the stress tensor) has to be equal to the mass carried by the 



using the identity \ — ^)'^ + A'^ = 6/2 and thereby obtain : 





(61) 
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particle number current Jj, i.e., 



Toi = mJi . 



(62) 



This is an example of algebraic identity between operators implementing sym- 
metries that hold in the microscopic theory and must be retained in the effec- 
tive theory. [13] 



Parametrization of the effective theory 

Throughout the analysis in this work, we shall use kf, (or, alternately, the 
spin-polarization P) and the pairing gap A as free input parameters that will 
be compared to experiments. As to this choice of free parameters, a remark 
is in order. In a truly microscopic theory, formulated in terms of the origi- 
nal fermions, the pairing gap would be determined by solving the gap equa- 
tion [43] together with equations specifying conservation of the total number 
of atoms and the population imbalance. The strength of coupling between 
fermions, naturally, shows up in these equations. Our theory, however, is not 
microscopic: we here assume the existence of the BPl phase and construct an 
effective field theory for this phase. Being formulated in terms of collective 
rather than microscopic degrees of freedom, our effective theory does not ex- 
plicitly have the coupling strength between original fermions and instead uses 
the pairing gap as an independent parameter. This choice is also motivated 
by the recent experimental developments in the field of atomic Fermi gases: 
it was demonstrated that using the rf-spectroscopy it is possible to measure 
the pairing gap by breaking fermion pairs. [15P6] Alternative methods of de- 
tecting a long-range pairing order in a degenerate Fermi gas have also been 
theoretically proposed, where the pairing function is directly measured in real 
space via a matter-wave interferometric techniques. j47j 




4-2 Symmetry-based effective Lagrangian 

The effective Lagrangian of the system ought to obey two global U{1) sym- 
metries, one of which corresponds to the total atom number conservation (to 
be denoted as Uc{l)), and the other one to the conservation of the difference 
in the number of atoms of spin- up and spin-down species (denoted as Us{l)). 
Our low-energy effective Lagrangian for the gapless branch of fermions (Bo- 
goliubov quasiparticles) , described by the field x(a;), and the superfiuid phase 
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field 6{x) (in the imaginary-time path-integral formalism, with r = it) reads 



C = x*[dr+eH^)]x+ci{dr9Y+C2{Vef+csX*X *5r^4 



1 



(V^)2j+V^.j+..., 
(63) 



2mp 



and represents an extension of the theory derived by Son and Stephanov [2D] 
to the case of an arbitrary spin population imbalance. In 0631) the ellipses 
stand for possible higher-order derivative terms of the 6 field; j = (x*Vx — 
Vx*x)/(2mpz) is the "paramagnetic" fermion (mass) current with rup = 2m 
being the total mass of the Cooper pair; e{—iV) is the operator form of the 
gapless fermion dispersion ([2]), written in the coordinate representation. The 
Lagrangian has the shift symmetry 6 — >■ 6+a, due to the Uc{l) particle number 
symmetry. Consequently, it contains the coordinate and time derivatives of 9, 
but not 9 itself. 

The phenomenological parameters ci, C2 and C3 are not constrained by the U{1) 
symmetries. While ci = dn/dfi {n being the total atomic density), C2 
and C3 are constrained by the superfiuid density n^. In this regard, an im- 
portant difference between the bosonic (phase) sector of our theory and the 
effective low-energy theories of bosonic superfiuids or neutral fully-gapped su- 
perconductors (with equal spin population) ought to be pointed out. Namely, 
in theories of the present type, in order to satisfy the Galilean invariance 
represented by the constraint (1^2]) . the low-energy effective Lagrangian can 
depend on the phase field only through the Galilean-invariant combination 
Ue = dr9 + 2^(^^)^ ("^0 being the mass of an elementary superfiuid con- 
stituent, e.g., the mass of a single atom in the case of '^He or the mass of a 
Cooper pair in case of fermionic superfiuids), that is. 



where P stands for an arbitrary polynomial. [l3] Keeping only the terms of the 
lowest order in the derivatives of 9, Ce reduces to the form ci{dr9Y + C2(V6')^, 
where coefficients ci and C2 are fixed by the requirement that this Lagrangian 
correctly describes the dynamics of the gapless Goldstone mode (Anderson- 
Bogohubov mode in the case of a neutral superconductor) associated with the 
spontaneously broken global U{1) symmetry. [Note that the term linear in 
dr9 is omitted, despite being of the lowest order, since it is a total derivative 
and the time-dependent topological configurations are not considered.] In our 
case, however, with additional low-energy degrees of freedom (gapless fermion 
excitations), the coefficient C2 is renormalized at every order of the effective 
theory and is constrained together with C3 by an additional requirement that 
the superfiuid density matches the one calculated from the microscopic theory. 
This identification will be made in the following section. 

The Galilean invariance of the fermion-dependent part of this Lagrangian is ex- 




(64) 
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plicitly demonstrated in Appendix |X1 An important consequence of this invari- 
ance is that the coefficient of the term V^-j must be unity. As a prerequisite for 
proving Gahlean invariance, we have shown that the Bogohubov-quasiparticle 
field remains invariant under Gahlean transformations. The transformation 
law for this quasiparticle field is thus essentially different from that of the 
original fermions, used as a basis for an alternative effective field theory of 
a polarized Fermi gas in Ref. [35]. This is consistent with a quite general 
argument that the transformation properties for quasiparticles in the low- 
energy effective theories should not depend on the quantities such as the 
bare particle mass m. [48j As a by-product of this transformation law, the 
Bogoliubov-quasiparticle current j is invariant under Galilean boosts, which 
is also consistent with the invariance of the quasiparticle momentum. 

4-3 Effective action for phase fluctuations 

Using the and a fields via Eq.(6), the Lagrangian (l63l) can be rewritten as 



In order to arrive at an effective phase-only action S[9] = S[(j), a], we integrate 
out the fermion field x '■ 



where S[x, 0] = S[x, <t>, a] = Jq dr J dii.C is the Euclidean action corresponding 
to Lagrangian (|65|) (with (3 = {ksT)'^ the inverse temperature). To this end, 
we first note that the fermion field enters Lagrangian (l65l) through a quadratic 
form x*Kx = X*i-Qo^ + X)x, where 



C = x*[dr + e{-iV)]x + ci (9.0)2 + C2(V0 - a)^ 
+csX*X ^dr<P + + (V0 - a) ■ j . 




(65) 




(66) 



1 



(67) 



is the noninteracting fermion propagator, and X = X'-^-* + X^^^ where 



X(^) = idrCj) + 



1 



(V0-a)- V , 



(68) 



1 



(V0-a)2 , 



(69) 



2 m, 



■p 



are respectively of the first and second orders in fields and a. 
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Integrating out the fermionic degrees of freedom gives rise to a contribution 
SF[(f),a] — — trlni^ to the effective action S[(f),a], where 

-trlnX = -trln(-e?o"^) - trln(l - QoX) . (70) 

The contribution of the self-energy X to the effective phase-only action is 
evaluated by employing the usual loop-expansion of the trace: by Taylor- 
expanding the second term on the right-hand side of the last equation (using 
ln(l — z) — — Z^^i z"'/n) we obtain 

OO 1 

- tr In ii- = const. + V - tr[(^o^)1 ■ (71) 
^1 ^ 

Using diagonality of the noninteracting fermion propagator in the momentum- 
frequency space {Qo{k,k') = Goik)5kk' with ^o(^) = (^'^n ~ ^^k)""*^, here dis- 
played using compact four- momentum notation : k = (k, ia;„)), it is straight- 
forward to show that 



tr(goX) = ^Y.^o(k)X,,, , (72) 
tr[(^?o^)'] = E Qo{k)Xk,k+Mk + q)Xk+,,k , (73) 

where Xk^k' stands for the Fourier transform of X. In order to obtain the 
effective action S[(j), a] to second order in fields 4> and a, we employ the above 
expansion to first order in X^'^^ (tree level) and to second order in X^^^ (one- 
loop order). 

The tree- level contribution of X^'^^ to S'^^f^.a] (and therefore to the effective 
phase-only action) is obtained by replacing x*X by its average value {x*x) ~ 
ni). It can easily be demonstrated that 

Xkl' = (^n - Un')(t)k-k' + TT— (k + k') • {afe_fe/ - i(k - k')0fe_fe/} , (74) 

where a^ = ak5uj„,o (the vortex gauge field is time-independent, i.e. classical). 
As a special case of the last equation, in the previously adopted Coulomb 
gauge (in which q • a±q = 0, hence q • a±g = 0) we obtain 



<-'=f'-i<»-(''+i)}^-'-;s;''-^-" P*" 

where q = (q, ia;;) and k_|_ = {(q x k) x q}/q^ is the transverse component of 
the three-dimensional vector k with respect to q. While it is easy to show that 
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the first order contribution (tree level) of X^^^ is equal to zero, by inserting 
the last two equations into Eq. fl75]) we find its contribution to S'i7'[(/), a] at 
one-loop order. 

The effective action for (f) and a is obtained by gathering S'i7'[(/),a] and the 
fermion-independent terms of the original action : 

a] = J^^ drj rfx [ciidr^f + C2(V0 - a)^] + a] . (77) 

In the momentum-frequency space, to second order in fields (p and a, it is 
represented by the quadratic form 



nb \ 2 1 / N / n(g)\ 2I , , 

C2 + 7^ C3 ) q + ^RiMmi + Ci - — - UJi 



q I 



2mp j 2m^ 



(78) 

The first term corresponds to the propagating Goldstone modes of broken 
U{1) symmetry, and the second one to its corresponding topological defects 
- vortices. [Summation over repeated indices in the last equation is implicit.] 
Here 

m = ^j:Qo{k)go{k + q) (79) 
is the fermion density polarization bubble, while 



RM = ^T.So{k)go{k + q){k,+ ^j){k, + ^f^ , (80) 

Piq) = ^T.So{k)go{k + q)^ (81) 

represent the longitudinal and transverse current-current correlation functions, 
respectively. In obtaining the form of the latter, we have made use of the 
identity 

EA:x.fcx,F(|k|) = ^EkW|k|) , (82) 
k k 

valid for any rotationally-invariant function F(|k|). The Matsubara frequency 

sum that is implicit in all of these response functions evaluates to 

T, 2^ GoikjGoik + q) = — — — . (83) 



It is important to point out that there is no RPA-type correction from the 
interaction vertex j • V0 to the transverse current-current correlation function; 
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this is manifest in our choice of the Coulomb gauge for the topological gauge 
field a. 

As proven in Appendix [Dl Rij{q) = R{q)6ij. Consequently, the phase-only 
action in Eq. (ITS]) in the zero-temperature static limit reduces to 



q \ 



where and are the zero-temperature static limits of P{q) and R{q), 
respectively. 

The superfiuid mass density ps, which plays the role of rigidity in the present 
problem ("spin- wave" stiffness in the XY-model terminology), |19] can be 
identified from the long-wavelength (q — > 0) limit through the relation 



This constraint on C2 and C3 can be equivalently stated as 



C2 + ^ C3 = ^ - ^q=" (86) 
2mp 2mp 



q=0 



and implies that the phase-only action in Eq. (!84|) adopts the form 



which is free of the phenomenological parameters of the original theory. 



5 Effective theory for vortices and the interaction potential 



Starting from the effective theory of phase fluctuations (described by (IHTI) ) 
and integrating out the regular (spin-wave) part of the phase field, we derive 
the effective action S'eff[a] for vortices : 

Along these lines, a straightforward Gaussian functional integration yields the 
result 

S-W-i:(£:; + ^^)a,-,. (89) 
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(Because the vortex gauge field belongs to the classical sector of the theory, 
the derived effective action contains only the ui = part). With the aid of 
identity (ITS]) , the last result can be conveniently recast as 



From the last equation we read off the momentum-space form of the effective 
interaction potential between the vortices : 



In addition to the long-range component proportional to 1/q^ (logarithmic 
interaction in the real space, i.e., 2D Coulomb potential), characteristic of 
the conventional two-dimensional charge-neutral superfluids, we have an ad- 
ditional component 



due to the presence of gapless fermions. 
5. 1 Properties of P^ and 

To calculate P^ one has to resort to a numerical evaluation. Yet, before em- 
barking on numerical work we can put P^ into a convenient analytical form. 
In Appendix [B] we demonstrate that P^ can be reduced to a two-dimensional 
principal- value integral 



p: 




b 



(93) 

where momenta k and q are expressed in units of kh and = |kp/2m — 
(fi/kf). The presence of the prefactor kf oc (n| — n^) indicates that in the 
thermodynamic limit the induced potential is proportional to the density of 
gapless fermions, which could have been expected on physical grounds. 

In the regime of small |q| (|q| < O.lk^) numerical evaluation becomes rather 
troublesome due to the strongly singular character of the integrand in Eq. (p3|) . 
However, as demonstrated in Appendix O) by replacing dispersion with its 
linearized form (et — Wb(|k| — kb)), for |q| ^ kb we can derive an expression 
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for Pq in the form of a controlled expansion in powers of |q|/A;fe : 



p2 



|q| 



Thus in the |q| limit we obtain: 



+ 




|q| < h 



(94) 



|q| 



(95) 



The last result can be given in a more concrete form. Applying the general 
expression 



dk 



k|=fct 



to the case of dispersion ([2]), we find 



Vb 



m 



A-2 



(96) 



(97) 



Inserting the last result into Eq. ( 195|) gives 



rnkf v \2m 



yU)2 + A2 



2m 



(|q| 



(9J 



Some typical results of numerical evaluation of the response function for 
0.1k/\ < |q| < /ca are displayed in Fig. ID. 41 (where 2kb < k^). The salient 
characteristic of these results is a knee- like feature at |q| = 2ki,, which reflects 
the existence of an effective Fermi surface with diameter 2kb. It bears analogy 
to the 2/i;ir-feature of the paramagnetic spin susceptibility in 3D, responsible 
for the RKKY indirect-exchange interaction between magnetic impurities in 
non-magnetic metals, [36] albeit the 2fcf,-feature found here comes from the 
current-current correlator so that it is different from the RKKY interaction in 
its dynamical origin. The values of obtained analytically in |q| — > limit, 
based on Eq. fl98l) . differ just slightly from numerical values at |q| = O.I/ca, 
indicating that P^ can be approximated as a constant in this numerically- 
inaccessible region < |q| < 0.1 A; a- The fact that P^ has very weak mo- 
mentum dependence at small q can be inferred from the coefficients in the 
controlled expansion of P^ given by Eq. flMl) . 

In Appendix [Dl using methodology analogous to the one employed in Ap- 
pendix O, we show that |q| — > limit of i?q is equal to that of P^ : 



pO 
q=0 



Gn'^Vb 



(99) 
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whereby Eq. (p2!) can be recast as 



2 pO _ pO 

V^ind(q) = ^^-^ . (100) 



Now, by virtue of controlled expansion (IMI) . we obtain that 

2m2 lOvr^tif, 



Wq) = -AT74r + ^(^) iH«h), (101) 



and, in particular. 









(|q| 



//•oo 
V^,d(r) dh = 271 rV,Ur) dr (102) 

is finite : 

l^nd(q = 0) = --^^-^ . (103) 
2m^ lOTT^ffo 



5.2 Effective vortex interaction potential in real space 



Let -F(|q|) be a rotationally- invariant function in momentum space and A the 
upper momentum cutoff. The inverse two-dimensional Fourier transform of 
F(|q|) is given by 

F(r) = -i-/ Fi\ci\) e^^-^ d\ . (104) 

(ZTTJ^ "'lql<A 

Using the identity 

e*l'il'^^°^'^c/y; = 27rJo(|q|r) , (105) 

Jo 

where Jo{x) is the zeroth-order Bessel function of the first kind, the last equa- 
tion becomes 

F{r) = ^ r Iqli'dqD'/odqIr) d\q\ . (106) 

ZTT JO 

In our effective theory, the upper momentum cutoff is set by A;a, thus the 
induced potential in real space is given by 

Mnd(r) = ^ |q|Mnd(|q|) Jo(|q|r) d\q\ , (107) 

ZTT JO 

viz., 

'^-W^S^i ^^MHr)M\. (108) 

Our numerical calculations of Mnd('") for different values of relevant parameters 
(/cb. A) show that the induced potential has damped oscillatory character. 
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closely resembling the spatial dependence of the RKKY exchange integral. As 
can be seen from Fig. ID. 51 this induced potential has alternating attractive 
{dVi-ciA/dr > 0) and repulsive {dVirid/dr < 0) parts. At short distances the 
induced potential is always attractive, and the first repulsive branch appears 
at the length scale r ~ (10 — 25) k^^, depending on the polarization. Spatial 
period of the observed oscillations is set by the spin polarization, but is not 
so simply related to the radius of the effective Fermi surface as in the case of 
genuine RKKY or Friedel oscillations. 

The total (effective) vortex-vortex interaction potential in real space is given 
by the sum of the induced potential and the conventional repulsive logarithmic 
potential. The latter is given by 



where the superfiuid density Ug is calculated in Sec. [31 As our calculations 
demonstrate, the effective vortex-vortex interaction shows three characteris- 
tic types of behavior, i.e. three polarization-dependent regimes. The critical 
polarizations corresponding to the boundaries between these different regimes 
are not universal but depend on the actual location in the part of the phase 
diagram pertaining to the BPl phase. 

In the regime of relatively low polarization, the total potential is dominated 
by the conventional repulsive logarithmic part; the effective vortex interaction 
is repulsive {dVes/dr < 0) at all distances. The resulting vortex phase is ac- 
cordingly expected to be conventional, with triangular vortex arrangement. 
An example is shown in Fig. ID. 61 

In the other extreme - the regime of high polarization, the induced potential 
plays a dominant role at short and intermediate distances. This renders the 
total potential attractive at short distances, with pronounced oscillating fea- 
tures resembling the RKKY interaction, as illustrated in Fig. ID. 71 Comparison 
of the induced and the total vortex-vortex interaction potential in the high 
polarization regime is depicted in Fig. ID. 81 

The attractive nature of two-body interaction already at short distances sug- 
gests an instability of the vortex lattice. However, whether this instability 
really occurs is still an open question for the following reasons. The physics 
at distances shorter than the healing length ^ (to be discussed in the next 
section) is not captured by our effective theory ; also, the multi-vortex in- 
teractions, not considered here but certainly allowed as higher orders in the 
effective vortex action, may support unusual vortex phases. This regime thus 
requires more elaborate further investigation. 

Apart from the two extreme regimes already described, in a narrow window 




(109) 
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of parameters the total potential is repulsive at short distances (r ^ (2 — 
3)A;^^) and becomes attractive at intermediate ones. This intermediate regime 
is illustrated in Fig. ID.9[ 



Due to the finite range of the RKKY-like induced potential, the truly long- 
distance dependence of the effective potential is governed by the infinite-range 
repulsive logarithmic interaction. However, for sufficiently large polarization, 
the effective potential is non-monotonous function of the distance between 
two vortices, a behavior that could potentially give rise to some exotic vortex- 
lattice structure. As is well known, the triangular-lattice configuration min- 
imizes the energy of a system of point vortices interacting through a repul- 
sive logarithmic (2D Coulomb) potential. [5U] Physically, this is a consequence 
of the fact that triangular vortex arrangement provides maximum nearest- 
neighbor distance at fixed vortex density per unit area, which is a natural 
tendency with purely repulsive interactions (at least in the continuum, i.e., 
in the absence of a vortex-pinning lattice structure). Examples of such be- 
havior can be found even in physical situations unrelated to vortices, such as 
the low-density limit of an electron gas, where a triangular Wigner crystal 
is formed. Interestingly, as can be inferred from Ref. [50], with the conven- 
tional logarithmic interaction the total energy of the triangular configuration 
of point-vortices is only around 0.8% smaller than that of the square-lattice 
configuration. Such a small difference, however, is not very surprising given 
that the lattice periods of these two configurations (for the same aerial vortex 



densities) are also not very different, namely atr = j2/^/3asq — 1.0746 a 




For our modified potential between vortices, which is not repulsive at all dis- 
tance scales, the structure of the vortex lattice is an open issue. It is worth 
mentioning that a non-monotonous interaction potential between vortices (al- 
beit without oscillating character) has recently been found in multicomponent 
superconductors by Babaev and Speight. [51] The authors have also predicted 
the existence of exotic (non-triangular) vortex-lattice structures. 



In general, the interactions between topological defects mediated by the en- 
vironment in which they are embedded is an important subject of current 
interest. Very interesting in this regard was the study of nodal-quasiparticle- 
induced interaction between vortices in d-wave superconductors performed by 
Nikolic and Sachdev. [52] They have found that the effect of quasiparticles 
on the effective vortex-vortex interaction (and, for that matter, some other 
properties of vortices) is not very dramatic. This can probably be ascribed 
to the nodal character of quasiparticle spectrum in d-wave superconductors, 
as compared to the fully-gapless situation that we are concerned with in the 
present work. 
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5.3 Experimental parameters and conditions 



In order to elucidate the realm of validity of our effective theory and make 
contact with experiments, it is useful to estimate the physical healing (coher- 
ence) length and compare it with the inverse of the momentum scale /ca- In 
this section, we discuss different regimes where our effective theory applies or 
may not be relevant. 

To that end, we analyze the bosonic sector of the theory. It is known from the 
BEC studies [53] that the healing length can be expressed as ^ = (Svma)"^/^, 
where n is bosonic density and a the corresponding scattering length. Generi- 
cally, this is the length scale set by the chemical potential of bosons {li^ / (2m^^) = 
^b)i expressed to lowest order in \fna?. Therefore, in our case the healing 
length can be expressed as ^ = (Svrnsam)"^''^, where is the molecular scat- 
tering length and the superfiuid density is the density of bosonic Feshbach 
molecules. Starting from expressions for ^ and k/\ we obtain 



^ V2mA J2meF{A/eF 



that is. 



/ A \ 1/2 

e 1 (A) (3vr^n)V3 

t-l /q^ / \ 1/2 / \ 1/2 



(110) 



(111) 



with ep = k\/{2m)^ where kp = (Stt^tt,)^/^ is the momentum scale set by the 
total fermion density. The last equation can be conveniently recast as 



1/2 

, (112) 



where 



:il3) 



is a dimensionless diluteness parameter. [20] [Recall the familiar results in 
three important limits: k, — oo (+oo) in the BEC (BCS) limit and k = at 
unitarity.] Eq. (11121) is equivalent to 



1/2 



0.6171 X ^ TP>\^\^'^ (114) 

and implies that ^ / A;^^ depends on three dimensionless ratios and the dilute- 
ness parameter. 
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To provide a quantitative estimate of the ratio ■C/^a ^^e parameter regime 
relevant for realization of the BPl state, it is useful to recall the relevant 
details of the mean-field phase-diagram of a polarized Fermi gas, based on the 
two-channel (i.e., boson- fermion) model. [T9|[2Ti| This mean-field theory yields 
quantitatively reliable results in the narrow-resonance regime, being exact in 
the limit of a vanishing resonance-width. [M] According to this phase diagram, 
for intermediate negative Feshbach-resonance detuning (z/) BPl (SFm) exists 
in the region between lines 6m ~ (boundary to unpolarized EEC superfluid) 
and 6ci ~ 1.3|i/| (the boundary to a phase separated state - the superfluid- 
normal coexistence region). [Note the following difference in notation: here, 
the chemical potential difference is denoted as 6, whereas in Ref. [T9lf55] it is 
2h while the detuning is denoted as 6.] 

It is known that in the strong-coupling BEC regime of a superfluid Fermi 
gas with equal populations of two hyperfine spin components (balanced Fermi 
mixtures) the molecular scattering length is given by am = 0.6 a/ (a/ being 
the scattering length between fermionic atoms). pS] For a polarized Fermi 
gas, however, as shown by Sheehy and Radzihovsky (see Eq. 16 and Fig. 
2 of Ref. |21j), the molecular scattering length decreases monotonously as 
a function of S/l^l (or, equivalently, of the polarization) and vanishes at the 
aforementioned boundary of first-order phase transition to the phase separated 
state. Therefore, as follows from Eq. (I114p . right at the boundary to phase 
separation and in the immediate vicinity of it the coherence length becomes 
much greater than k'^^, thus making the quantitative implications of our theory 
not directly applicable in this special case. 

Taking am{S = 0) in place of am, together with typical values of Us/n and A/ep 
{ep = kp/{2m), where kp = (Svr^n)^/^ is the momentum scale set by the total 
fermion density) in the BEC regime, we estimate that ^ is of the same order 
as fc^^ when \k\ ~ 1-100. From numerical results for am{\j^\,S)/am{\j^\,S = 
0) obtained in Ref. [21], on account of the fact that ^/k'^ oc {am/ cif)~^^'^ ■, 
we can infer that the above estimate is just slightly modified as a result of 
am decreasing as a function of for example, for 5/\v\ = 1.2 the true 

molecular scattering length is an order of magnitude smaller than that of 
the unpolarized system, but the ratio i/k'^ is modified only by a factor of 
a/To ~ 3.16. For smaller values of 5 this factor is even smaller, i.e., it is of 
the order of unity. Thus, this estimate confirms that our choice of k/^ as the 
upper momentum cutoff of the theory is physically pertinent. 

Moreover, using the expression of A/e^ in Ref. 0, we can straightforwardly 
infer that in the BEC limit ^/k'^ oc IkI-"^/^. It follows that ^/k^ — > oo in the 
BEC limit, which seems to suggest that this limit is out of the application 
scope of our theory, since the latter is intrinsically valid for physics at dis- 
tances longer than ^. It is, however, important to emphasize that the effect 
of gapless fermions on the interaction between vortices is not even expected 
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to bear any physical relevance in the BEC limit, where the system at hand 
essentially becomes a Bose-Fermi mixture akin to the ^He-'^He mixture. [35] 
Technically speaking, this point is manifest in the Nishida-Son formulation of 
the effective Lagrangian for the imbalanced Fermi gas, through a vanishing 
coupling coefficient between the fermion current and the gradient of phase field 
(super current) in this limit. In our case, largeness of the physically-allowed 
inter-vortex distance scale (r > ^) compared to /c^^ in the BEC limit and the 
fact that at very long distances (compared to fc^^) the fermion-induced part of 
the vortex-vortex interaction is quantitatively unimportant compared to the 
conventional (infinitely-ranged) repulsive logarithmic contribution are indeed 
suggesting that the physical effect under consideration is absent in this limit. 
This is an important consistency check of our results. 



6 Summary and conclusions 

In summary, starting from a Lagrangian for the superfiuid phase field and the 
gapless branch of fermionic quasiparticles, we have obtained the effective ac- 
tion for vortices in a spin-polarized homogeneous superfiuid state with a single 
gapless Fermi surface. We have demonstrated that besides the conventional 
repulsive logarithmic part (2D Coulomb potential) the effective vortex inter- 
action potential has an additional, predominantly attractive, component in- 
duced by the presence of gapless fermions. This fermion-induced potential has 
oscillating character analogous to the RKKY indirect-exchange interaction. 
Interactions between defects mediated by the continuum they are immersed 
in (either bosonic or fermionic) have been studied quite recently in several 
different physical contexts and different dimensionalities. [FTIISS] Our work, 
however, constitutes the first study of this kind that concerns the interaction 
between vortices in superfiuids. It shows that besides the Friedel oscillations 
(charge sector) and the RKKY (spin sector), an analogous oscillating phe- 
nomenon appears in the vortex sector. 

Our study opens up a question as to the nature of the vortex lattice in gap- 
less fermionic superfiuids. Due to the partly attractive nature of the effective 
vortex potential that we have found, the resulting vortex lattice structure 
in BPl superfiuid phase could be different than the triangular lattice, which 
would be a spectacular experimental signature. The complexity of the problem, 
however, calls for an elaborate future study. Even when the potential has a 
unique distance dependence of a known analytical form (such as, for example, 
the conventional logarithmic interaction), calculation of the resulting lattice 
structure is quite a nontrivial task, since lattice-summation-methods [501159] 
for long-range potentials are strongly dependent on the actual form of the po- 
tential. The new vortex- vortex interaction potential is not obtained, due to the 
complexity of the problem, in a closed analytical form and has both attractive 
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and repulsive parts. This unusual, non-monotonous, distance dependence of 
the effective potential implies that vortex lattice structure may in fact not be 
unique, but also depend on the geometrical constraints on the system, for ex- 
ample, the range of distances between individual vortices realized for a given 
size of the superfluid container. The standard lattice summation methods may 
not be applicable and more sophisticated strategies need to be employed, for 
instance Monte Carlo calculations. 

In the present work we have studied the intrinsic effect of gapless fermionic 
excitations on the interaction between vortices in the BPl state and have 
therefore considered only the homogeneous case. Our results are expected 
to be also valid for a trapped system as long as the trap potential varies 
smoothly on the scale of the Fermi wavelength (or, more generally, the longest 
physical length-scale in the problem), that is, in the regime of validity of 
the local density approximation. However, an important problem yet to be 
explored is the possible influence of strong spatial inhomogeneities caused by 
the presence of the trap on the form of the vortex lattice, as studied by Sheehy 
and Radzihovsky in the context of trapped Bose gases. |60,61| They have 
provided an explanation for the striking uniformity of the vortex lattices seen 
in experiments in spite of the strong spatial variation of the local superfluid 
density imposed by the trap. Moreover, they have shown that an interplay of an 
inhomogeneous trap potential and vortex discreteness leads to a vortex density 
that is largest in the center of the trap, a counterintuitive result from the 
energetic point of view because both the kinetic energy cost and the repulsive 
interaction between vortices are proportional to the local superfluid density 
and are therefore largest in the center of the trap. As we have shown in the 
present study of a spin-polarized Fermi gas, for sufficiently high polarization 
the effective interaction between vortices in this system is attractive at short 
distances and could therefore bring about some completely new effects, such as 
the competition between this attractive interaction and the kinetic energy cost. 
Further investigation of the properties of "vortex matter" in spin-polarized 
Fermi gases is thus clearly called for. 
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A Galilean invariance 



In this Appendix, we explicitly demonstrate the Galilean invariance of the 
fermion-dependent part of Lagrangian fl63|) . We shall first establish an explicit 
relation of the quasiparticle field in the present effective field theory to the 
fermion particle field in a microscopic model, and then derive the Galilean 
transformation properties of the quasiparticle field from that of the (micro- 
scopic) fermion fields. Subsequently, an alternative approach will be given to 
provide a further justification and understanding. 

A.l Microscopic relation of the quasiparticle field 

To examine the Galilean transformation of the Bogoliubov quasiparticle field, 
let us consider a microscopic model of Lagrangian (fT9|) . As a result of the 
Hubbard-Stratonovich transformation in the Cooper channel, introducing the 
auxiliary pair field A(x), this Lagrangian changes to 

>C = - ^ - /i.) + (^T^IA(x) + c.c) + -^\A{x)\' . (A.l) 

[Summation over repeated pseudo-spin indices in the last equation is implicit.] 
Ignoring fluctuations of the amplitude of the order parameter, i.e., assuming 
that A(x) = Ae^^^^\ it is advantageous to transform the fermion fields at each 
space-time point as [33] 

M^) = Mx) e^^(^) , rA^) = C(x) e-^^(^) . (A.2) 

This local (gauge) transformation is designed to transform away the phase- 
fluctuation dependence from the off-diagonal pairing potential terms to the 
diagonal (kinematic) terms in the fermion sector of the theory. As a result, 
the ipcj fermion fields are locally stripped off of any dependence on the U(l) 
phase 9{x). The transformed Lagrangian can be written as C = Cq + C^q, 
where 

is the mean-field Lagrangian for ipa- fermions, and 

^^,e = rja (tdrO + ^(V^^)') - ^ (CVV^. - VrJa) " . (A.4) 

With Cq naturally giving rise to the Bogoliubov quasiparticles as its elemen- 
tary excitations, q essentially contains, in an implicit form, all the couplings 
of these excitations to the superfluid phase fluctuations. 
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A cautious remark is needed for the Lagrangian derived above. It appears 
that we have just provided a derivation for the postulated effective Lagrangian 
(1631) . One may be tempted to determine the "phenomenological" coefficients 
of this Lagrangian in this way. For weak couphng, this can indeed be done. For 
a strongly interacting Fermi gas, the derivation from the microscopic model 
cannot be done in a controlled approximation, once the pairing amplitude 
and density fluctuations are included. The symmetry-based Lagrangian of the 
postulated form fl63l) describes the same physics, albeit from a more phe- 
nomenological point of view. Moreover, it does not suffer from the difficulty 
in strong coupling. In summary, the above derivation is understood to provide 
an example of how to separate the low energy Goldstone bosons (the phase 
fluctuation) from other degrees of freedom, but not a rigorous proof of the 
effective Lagrangian itself on a microscopic level. 

A. 2 Galilean transformation for quasiparticles 

Let us denote the laboratory frame as K and the corresponding spatial and 
time coordinates as x and t. We shall also denote a frame moving with velocity 
u relative to K as K', and its spatial and time coordinates as x' and t'. Un- 
der Galilean boost transformation with velocity u the space-time coordinates 
transform as : x — > x' = x — ut, t ^ t' = t. The spatial and time derivatives 
transform as 

V — . V' = V , dr — >d'^, =dr-i{u-V) . (A.5) 

Under the above Galilean transformation, the (microscopic) fermion field of 
mass m transforms in the standard way, 

ij{x) — > tp'ix') = e*(-™+5™"'*)^(a;) (A.6) 

where a; = (x, r = it). Being locally stripped off of any 'charge' Uc(l) phase 
dependence, the ipa fermion fields are by construction invariant under Galilean 
transformation: 

Mx) ^U^') = i>M , (A.7) 
and the 'charge' Uc(l) phase transforms 

e{x) — > e'{x') = e{x) - 2mu ■ x - zmuV . (A.8) 

The Bogoliubov quasiparticles (the two branches being denoted by xt Xi) 
can now be introduced as 

XT(k,r)=Mk^l(k,r)+t;k^J(-k,r) , (A.9) 
Xi(k,r)=M_kVii(k,r) -t;_k^]:(-k,r) , (A.IO) 



29 



where Wk, f k are Bogohubov amphtudes. For the lower branch quasiparticle 
field, which is simply denoted by x (i-^, X = Xi): this becomes in real space 



(A.ll) 



Xlx, T)= I dy [u{y - x)^/'|(y, r) - w(x - y)^/j*{y, r) 
By using (]A.7p . we find 

X{x) — > x{x') = x{x) . 



(A.12) 



From the transformation properties of x and 6, it is straightforward to prove 
the Galilean-invariance of the term 



X*X \ idr9 + 



1 



2m, 



X*xUe 



(A.13) 



in Lagrangian ( l63l) . Besides, using relations (1A.8I1 and ( lA.Sj) (the latter implies 
that e{—iV) is an invariant), we can readily prove that the combination 



x*[dr + e{-t\/)]x + ve-i 



(A. 14) 



is Galilean invariant up to an unimportant total derivative. This proves the 
Galilean invariance for the fermion- dependent part of Lagrangian (!63l) . 



A. 3 The Doppler shift 



An alternative check of the Galilean invariance of the quasiparticle field can 
be obtained by starting from the requirement that the quasiparticle energy is 
Doppler shifted under a Galilean boost. 

We first review a standard derivation of the Galilean transformation |62j . 
Recall how momentum and energy of particles with quadratic dispersion (e.g., 
bare fermions) transform under this Galilean boost : 

p — > p' = p - mu , E — > E' = E - p-u + ^mu^ . (A.15) 

Using these rules, it is straightforward to show that the combination p ■ x — E't 
shifts by a factor — mu ■ x + ^mu^t, which depends on the parameters of the 
transformation (m, u) but does not depend on p. Accordingly, every plane 
wave 

V7p (x, t) = const X e'(P-"-^*) (A. 16) 

acquires the same phase factor exp[i(— mu ■ x + |mu^t)] under the Galilean 
boost, regardless of p. Moreover, since an arbitrary single particle wave- 
function can be expanded in plane- waves (lA.16p . we conclude that each wave 
function picks up that same phase factor under this boost. Because Galilean 
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transformations are space-time symmetry transformations, the transformation 
property of the single-particle wave-function carries over to the field operators 

^(x,t) = ^a„0„(x,t) , (A.17) 

n 

where (f)n{'^,t) form an arbitrary complete orthonormal set of single particle 
states. 

Bearing in mind the definition ([S]), as a by-product of the transformation rule 
found above, we conclude that the superfluid phase field is transformed as ^62] 

^(x, t) — > e\x', t') = 9{x, t) - nipU ■ X + ^mpU^t , (A. 18) 

that is, we recover transformation (1A.8P when changing over to the imaginary 
time. By making use of the transformation properties lA.Sl it is straightforward 
to show that the combination 

Ug = id^e + ^{vef (A. 19) 

2mp 

is invariant under the Galilean transformation. 

In order to determine how the Bogoliubov quasiparticle field transforms under 
the Galilean transformation, we recall that the momentum of a quasiparticle 
is invariant under the Galilean transformation while the quasiparticle energy 
is Doppler-shifted (to leading order in the boost velocity): |3S] 

p — >p' = p, E — > E' = E -p-u . (A.20) 

Based on these properties, it is easy to demonstrate that for Bogoliubov quasi- 
particles the combination p ■ x — Et remains invariant under the Galilean 
transformations (independent of p and E), which using analogous reason- 
ing as above implies that an arbitrary single-particle wave-function and the 
field operator x(x;^) of ^ Bogoliubov quasiparticle is invariant under Galilean 
transformations : 

x(x,t)^X'(x',t') = x(x,t) . (A.21) 
This is the equivalent form of (IA.12P in operator formalism. 



B Expression for the transverse current response function 

Because > for |k| > kb and np{e) — * 0{—e) as T ^ 0, in the zero- 
temperature static limit the response function P{q) (defined by Eq. (15T1) ) 
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reduces to 



(27r)3 



^^(|k + q| - h)9{h - |k|) 9{h - |k + q\)9{\k\ - h 



IT] 



^(B.l) 

where t] 0+ and the momentum sum in (IHTj) has been replaced by an 
integraL We now undertake the change of variables k' = — k — q in the second 
term of the last equation. Because k'_|_ = — k_|_, we have that k'^ = k^. In 
other words, k^^ is invariant under this change of variables. Consequently, we 
arrive at 



(2^ 



|k + q|-fc,)^(A;,-|k| 



1 



1 



B.2) 

where the superfluous prime has been omitted (i.e., we have returned to the 
initial integration variable k). The last equation can obviously be simplified 

to 



P° 



^(|k + q|-A;5)^^(A;,-|k|) 



By virtue of the Sohotsky-Plemelj formula 



;^^0+) . (B.3) 



lim 

v\o X ±irj 



X 



(B.4) 



we can now demonstrate that Im {P°} = and that Re {P°} = P° IS given 



by 



V 



d^k 



|k + q|-A;5)^(A;6-|k|) 



(27r)^ ■■ - £:k - ^k+q 
where V stands for the Cauchy principal value. Using the identity d{x 
1 — 6{—x) for a; = |k + q| — /c;,, the last equation becomes 



(B.5) 



P" = P 



d^k 
(2^ 



l-^^(A;,-|k + q|) 



e{h-\k\) 



ki 



£^k ~ ^k+q 



(B.6) 



The term that contains the product of two step functions vanishes identically 
after the integration, because this product is even under the interchange k ^ 
k + q, while the fraction k^/ (ffk — ^k+q) is odd under the same transformation 
(while £k ~ ^k+q is obviously odd, the fact that k^ = (k + q)^ implies that 
ki is even) ; accordingly, we have 



V 



d^k 

(27r)3 



^(A:6-|k|) 



£^k — £^kH 



(B.7) 



With the aid of identity k^ = |kp(l— cos^ 9) and momentum re-scaling k/kh — 
k (such that all momenta are expressed in units of k}y) we express P^ as a 
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principal- value integral over the dimensionless momentum : 

K = jiv2'^ n^\'d\\-\ r'-^^sinede . (B.8) 

(27r)^ Jo Jo - ^k+q 

Finally, upon inserting dispersion ([2]) and making substitution x = cos 6, this 
integral leads to 

= 7%'^ / Ikr^lkl dx- 



(27r)2 ./n ' ' ' './-I / TT^ r 



-1 



(B.9) 



where = |kp/2m — {ii/kl) 



C Behavior of for |q| ^ fcf, and the |q| — > limit 



The most general expression for P{q) = P{(\,iuJi) reads 

^(q, ^^i) = ^ E —7—7 — k^ ' 



that is 



P(q, luj,) = - / 7^,^:^ kl . (C.2) 

2 7 (27r)^ lui + Sk - £k+q 



In order to calculate P^ for |q| ^ fcfe, we start from the expansion 



n^(^k) - ri^(5k+q) = ^^^^ (£k - ek+q) + 0{\q\') , (C.3) 

C^k 

valid for |q| ^ /c;,. At zero temperature ^^^(e) = 6{~6), implying that dnp{e)/de 
—S{e). For linearized dispersion = f^bilM — kb), using the fact that 6{cx) = 
6{x)/\c\, we find 

npie^) - riF(£k+q) = (|k + q| - |k|) 5{\k\ - h) + OdqH . (C.4) 

Here k ■ q = |k||q|cos6', and consequently |k + q| = (|kp + 2|k||q|cos6' + 
|qnV2. 



On account of result flC.4p . together with k^ = |kp(l — cos^ 9), Eq. flC.2l) leads 



to an integral (trivial integration over the azimuthal angle yields factor 27r) 

1 roo rl 

J- / 1/1 , / 



I /■oo ri 



(l-cos2^)(|k + q 




kmk\-h) 


iu;, - t;fe(|k + q 




k|) 



(C.5) 
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Upon executing the integral over |k| and introducing substitution x = cos 9, 
we arrive at 







q X + 


q 




iuji - Vb 




q X + 


q 





Piq,tcoi) ~ — ^ / ;\ f dx . (C.6) 



Another variable substitution t = yk^ + 2kb\q\x + |qP turns the last integral 
into 

P(q, ~ / — dt . 

(C.7) 

By carrying out this integral and taking the static limit Ui —>■ 0, we obtain the 
result (without the prefactor) — yA;^|qp, implying that the first order term in 
expansion (IC.SP yields the q-independent contribution 

k^ 

OTT'^Vb 

to Pq. In a similar manner, lengthy but otherwise straightforward calculation 
shows that the next (second-order) term in expansion (IC.Sp . namely 



(irt-ek+,)' = -r(|k + q|-|k|)V(|k|-fc) , (C.9) 



'k 

adds the contribution 

ki , ,0 1 



q| +7T^7r^|qr ■ (cio) 



Therefore, for |q| -C fc^ this response function is given by 



6 

implying that 



^ Gn^Vb lOn^Vb 



(|q|-o). (C.12) 



D Calculation of -R° (q) in the |q| ^ limit 



The most general expression for Rij{q) = Rij{q,iLJi) reads 
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that is 



' (27r)^ za;i + ek-£k+q V 2 J \ ' 2' 



We first show that Rij{c[, iuji) = for i ^ j. To that end, we perform a rotation 
of the coordinate system around the z-axis that maps the x-axis onto the y- 
axis and the y-axis onto -x. Knowing that the module of the jacobian of this 
transformation (rotation) is unity and that depends only on |k| (which is 
invariant under this transformation) we obtain that Rxy{q) = —Rxyil) and 
Ryx{q) = -Ryx{q), which implies that Rxyiq) = Ryx{q) = 0. 

In order to calculate 

we perform a rotation of the coordinate system that maps the z-axis onto the 
z-axis, while leaving the remaining axis invariant. i?jj(q, zw;) then becomes 

J (27r)^ tuji + e-k' - e-k'+d \ 2J 

for both i = X and i = y. Thus Rxx{q) = Ryy{q) = R{q), and since the last 
integral can depend only on |q| we can choose q to lie along the z-axis, in 
which case R{q) can be expressed as 

R{q, lui) = / . — + — , p.5 

J (27r)^ icui + £k' - ^k'+q V 2y 

i.e. as 

R{q,iuji) = ——— ■ ^ kcos0 + — . p.6 

J (27r)3 tUJi+Ek-Ek+d \ 2 



By employing transformations analogous to flC.3l) - flC.7p in the calculation of 
Pq we arrive at the expression for R{q,iuJi) in the zero-temperature limit : 

Rlq.iuJi) c::^ — — — — — -r / ; — dt . (D-7) 

8(27r)2t;,|q|3 t - h - ^ ^ 



By carrying out this integral and taking the static limit cj; — 0, we obtain 



K^-Ar (lql-o). (D.8) 
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FIG. D.2. Superfluid density in units of the total atomic density n = n^+n-^ for three 
different values of the chemical potential mismatch 6. Both 6 and A are expressed 
in units of |/i|. 
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FIG. D.5. Induced vortex interaction potential in real space (in units of A), for 
m = 1.0 and A/|/i| = 1.0. Values of kf, are given in units of k/\: kf, = 1.352; 1.640 
correspond to polarizations P = 0.702; 0.808, respectively. 
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FIG. D.6. Effective vortex interaction potential in real space (in units of A), for 
m = 1.0 and A/|/i| = 2.0. Values of kf, are given in units of k/^: kj, = 0.623; 0.724; 
0.826 correspond to polarizations P = 0.155; 0.227; 0.314, respectively. 
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FIG. D.7. Effective vortex interaction potential in real space (in units of A), for 
m = 1.0 and A/|/x| = 1.0. Values of A;^ are given in units of k^: kj, = 1.500; 1.352; 
1.205 correspond to polarizations P = 0.763; 0.702; 0.626, respectively. 
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FIG. D.8. Comparison of the induced and the effective vortex interaction potential 
in real space (in units of A), for m = 1.0, A/|/i| = 2.0. Values of kh (in units of /ca) 
and P are indicated in the plot. 
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FIG. D.9. Effective vortex interaction potential in real space (in units of A), for 
m = 1.0 and A/|/x| = 1.0. kf, = 1.022 (in units of A;a) corresponds to P = 0.504. 
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